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1 Dual systems 

Study of the nonlinear lattices 

Qn = r(q n )(g n+ i - g n ), g n = g(q n - Qn-i) (!) 

have brought Toda to the notion of dual system, which is obtained from the 
original one by "replacement of particles by springs and springs by particles 
accordingly to the certain rules" ij. These rules are the following. The lattice 
(0) is equivalent to the Euler-Lagrange equations for the functional 



C = J dtJ2K<in)-K), a" = l/r, b' 



9, 



where a and b define kinetic and potential energies respectively. The momentum 
conservation law (T n denotes shift operator q n — > q n +i) 

^a'(<7„) = (T n - l)g n 
allows to introduce the new variable Q n accordingly to the formulae 

Qn = gn, Qn+l — Qn = a '(<?ra)- 

Following the terminology of the paper |2j , we call this and analogous changes 
generalized Legendre transformations. The variable Q n satisfies the dual lattice 

Qn = R{Qn){G n +l — Gn) 
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where the functions R and G are defined by equations g' — R{g), G(a'(q)) = q. 
For example, the Toda lattice q n = exp(g n +i — q n ) — exp(g n — q n -i) is dual to 
the lattice Q n — Q n (Qn+i — 2Q„ - Q n -x)- 

The following features of the considered transformation should be mentioned. 
Firstly, it is possible due to the invariance of the Lagrangian with respect to 
the shift q —* q + const . Secondly, it can be performed for arbitrary r and g, 
and therefore it is not related with integrability. At last, it is involutory (up to 
the shift T„) and therefore it cannot be used for reproducing of solutions of the 
lattice. 

It was shown in B that the situation changes drastically in the case of the 
relativistic Toda type lattices 

'An = r(q n )(q n +ih n+1 - q n -ih n + g n +i ~ 9n)- (2) 

The Lagrangian for the lattice (||) is obtained from the one written above by 
adding of the term —q n c n , d = h which is interpreted as the presence of magnetic 
field. This term is equivalent, up to the total divergence, to —q n c n +i, and this 
leads to two equivalent forms of the momentum conservation law 

^-(a'(q n ) - c n ) = (T n - l)(q n h n + g n ) ^> 
at 

^( a '(9«) _ c «+i) = i T n - l){q n -ih n + g n ) 
and consequently to the pair of generalized Legendre transformations 

T— . Qn qnh>n ~t~ Qn j Qn+\ Qn — Q Cn j 

T+ '■ Q n = Qn-lh-n + 9m Qn+l — Qn = a {<ln) — c n+l- 

In contrast to the previous example, these transformations generally lead out 
off the class (|^). The requirement that the transformed lattice must be of the 
same form turned out to be severe enough to isolate exactly the subclass of 
the integrable lattices (that is, possessing higher symmetries and conservation 
laws) 0j. It is not so surprisingly, because in this case the combination TZ 1 T + 
of the Legendre transformations defines the Backlund autotransformation of 
the lattice under consideration, and this property is indispensable feature of 
any integrable system. The situation is quite analogous to the case of the 
Backlund transformation for the KdV equation which is composition of two 
slightly different Miura transformations. 

The aim of the present paper is to apply these ideas in the totally discrete 
case: the differentiation d/dt is replaced by the shift T m : q rn ^ n — * Qm+i,n on 
the second discrete variable. In order to abridge notation we will omit, as a 
rule, the subscripts m,n : q — q m . n , qi,o — q m +i.n and so on. Moreover, it is 
convenient to introduce notations for the differences (see figure below) 

x = q-q-i,a, y = q-q -i, z = q-q- X - X . 
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Obviously the following identity is valid 

(T m - l)yo,i = (r„ - l)xi i0 . (3) 
The analog of the lattice (Q) is given by 5-point difference equation 

(T m - l)f(x) = (T n - l)g(y) (4) 

of the discrete Toda lattice type || [|, || , and Legendre transformation is of 
the form 

*o,i = Q0.1 - Q = f{x), X lf0 = Qx,q - Q = g(y). 

It is not related to integrability as well as in continuous case. The analog of the 
relativistic lattices (|^) is more symmetric than its continuous counterpart. Let 
us consider the functionals of the form 

C = + b(y) + c(z)). 

Obviously, the rights of all terms (which play the roles of the kinetic and poten- 
tial energies and magnetic field) are absolutely equal, and therefore it is more 
aesthetically to enumerate the variables q by the points of the rectilinear tri- 
angular lattice rather than square one. Nevertheless we use the standard grid 
which makes computations more easy. Euler equation is of the form 

(T m - l)f{x) + (T„ - 1)5(2/) + (T m T„ - l)h(z) = (5) 

where / = a' , g — b' , h — c' and we assume that f'g'h' =/= in order to exclude 
the case (^). Up to the author knowledge, the equations of this type appeared 
for the first time in the papers || [?], ||] by Suris, who also obtained them starting 
from the relativistic Toda lattices. 
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The third difference can be expressed through x, y in two different ways and 
this circumstance allows to rewrite (||) in two equivalent forms of the momentum 
conservation law: 

(T m - l)(f(x) + h(z 0A )) + (T„ - l)(g(y) + h{z)) = & 

(T m - l)(/(a;) + h{z)) + (T„ - l)( 5 (y) + /i(z li0 )) = 0. 
As a result we can define a pair of the Legendre transformations 

T+: (X,Y_ 1>0 ) = T(x, 2/ 0) i), T_ : (X ,-i, F) = Tfri.o, tf) (6) 

where the mapping T : (x, y) — > (X, y) is given by formulae 

X = g(y) + h(x + y), Y = —f(x) — h(x + y). (7) 

Notice, that the lattices (||) are the continuous limit of the equations (||). 
Indeed, let us consider the family of the functionals 

£ £ = y^(ga(a: m ,n/£) + eb(y m . n ) - c{y m .n) + c(z m ,n)), 

m,n 

and let q m ^ n = q n (t), t = me. It is easy to see that passage to the limit e — > 
brings to the functional 



which coincide with the one considered above, up to designations. Legendre 
transformations for the corresponding relativistic Toda lattice also are obtained 
from @ by the mentioned limit. 

Definition 1 Equation (^) is called integrable if the generalized Legendre trans- 
formations (Q) are invertible ( on the differences ) and map it into the equa- 
tion of the same type. 

Theorem 1 Equation ^) is integrable in the sense of the Definition 1 iff the 
inverse of (j^j is of the form 

x = G(Y) + H(X + Y), y = -F(X)-H(X + Y). (8) 

In this case the dual equation is 

(T m - 1)F(X) + (T„ - l)G(F) + (T m T n - l)H(Z) = 0. (9) 
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Proof. Obviously, it is sufficient to consider only transformation T+ . Let T^T 
be of the form stated above, then one can easily check that identity (^) is 
equivalent to the equation (||). Conversely, assume that T + maps (||) into ([^). 
Solving the formulae (X, i^-i.o) — T(x,yo,i) with respect to x,yo,i one obtains 

x = $(X,Y_ 1<0 ), yo.i = *(*,y_i, ) 
and identity (Q) yields equation 

$(^i,i, K ,i) - $(*i,o, y) - *(x 1)0 , y) + y_i, ) - o 

which must be equivalent to (^|), that is 

F(Xi, ) - f(x) + G(y ,i) - G(y) + H{x 1A + y 0; i) - + y_i, ) = o. 

Let us consider Xi t ± as the function on the rest variables involved in these 
equations, then 

dX hl _ g$ , d<5> _ G'(Y . 1 ) + H'jXj.j +y ,x) 
Sy ,i ~ 8Y ,i 1 dX 1A ~ H'(X hl + Y 0t i) 

and hence $(X, y) = <^(G(y) + H(X + Y)). Analogously the equality 

ay-i, _ a* _ + g'pr + y_i, ) 
~ dX' 8Y- lfi ~ H'(X + y_i, ) 

yields y) = ip(F(X) + + y)). Finally, the relation 

~ dY BY I dX lfi ~ F'(X lfi ) 

proves Y) + &(X, Y) = \{F(X) - G(Y)). Using these relations one easily 
obtains ip(z) = az+{3, ip(z) — —az+j and finish the proof (notice that functions 
F, G, H from @ are defined up to the linear transformation) . I 

The immediate corollary is that equation which is dual to the dual equation 
coincides with the original one. This allows to use composition TZ X T + of the 
Legendre transformation for reproducing of the solutions. 

2 Classification theorem 

So, in order to classify the integrable equations (||) it is sufficient to find all 
functions f,g,h such that inverse of the transformation (^) is given by (||). 
Hence the Jacobian A = f'g' + g'h! + h'f of the map (0) must be nonzero and 
the following identities must hold 

XXY + VXY = 0, X X Y = X X X, VXY = VYY ■ 
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These three conditions are equivalent. Indeed, the Jacobi matrix is 

-h! f + hi 
-g' - h' b! 

that is xx = —yy- Straightforward computation proves that functions f(x), 
g(y), h{x + y) must satisfy the equation 

f" a" h" 

(9' + tiyj + (f' + h') 9 7 = (f+g>)^. 

Designations /' = 1/u, g' = l/v, h! = l/w rewrite it in more convenient form 

[v(y) + w(x + y)]u'(x) + [u(x) + w(x + y)]v'(y) = [u(x) + v(y)]w'(x + y). (10) 

The classification problem is reduced to solving of this functional equation. Of 
course, the functions u, v, w can be multiplied on the arbitrary constant and 
linear transformation 

x = c(x - x ), y = c(y - y ) <t=S> q m ^ n = c(q„ hn - mx - ny - const) 

can be applied, as well as permutation of the coordinate axes. 

At first let us consider some degenerate cases. Assume that two of the 
functions u,v,w are constant, say u and v. Then ( |To| ) yields that either w is 
constant as well or u = — v and w is arbitrary. In the first case the equation (^|) 
is linear and in the second it is of the form 

«(<?i,o + 5-1,0 - 5o,i - 9o, -1) + M<7i,i ~ 0) ~ HQ - 9-1, -l) = 
and admits one integration 

Ot(q m +l,n — Qm,n+l) + ^(ffm+l.n+l — Qn 



lni.il I 



Other degenerate case corresponds to the case of zero Jacobian, that is u + v + 
w = (this implies that all three functions arc linear). In this case one can prove 
that equation (^) can be reduced to the equation on the variable p = yo,i/x. 
Further on we will not consider these degenerate cases. 

Lemma 1 Functions u, v, w satisfy equations 

{u'f = 8u 2 + 2au + e, {v'f = Sv 2 + 2(3v + e, {w'f = Sw 2 + 2 1W + e. (11) 

Proof. At first prove that functions u{x) and v(y) satisfy equation 

(u" -v")(u + v)-(u') 2 + (v') 2 = k(u-v), k = const. (12) 



G 



Let us eliminate w from (|l0|). Applying the operator d x — d y one obtains the 
linear system on w, w' : 



w + it 



(u — v) 



Its determinant A is exactly the left hand side of the equation (|12|). If it is 
identically zero then ( p"2| ) is proved, otherwise one finds 

I U — V u , 2 /, , ,x ; U ~ V , II I . I II \ 

w + u = — — — (uu — (u ) + U V + U V ), w — — — — (u V + U V ) 

and consequently ((u — v) / A) y {uu" — (u 1 ) 2 + u"v + u'v') = 0. Assume that 
the expression in the second bracket vanishes. If u' ^ then v' — —u"v/u' + 
u' — uu" /u' , v" = —u"v'/u\ but then, as one easily checks, A = 0. If u' = Q 
then w + u = that is we come to the degenerate solution excluded above. 
Therefore ((it — v)/A) y = 0. Due to the symmetry between u and v one can 
prove analogously ((it — v)/A) x = and obtain (p^). 
Further on, rewriting (|12| ) in the form 

A2 



v" + k {v'f + 2kv 



u-\-vJ u + v (u + v) 2 
multiplying by v! /(u + v) and integrating with respect to x yield 
(u') 2 = 5(y)(u + v) 2 - 2(v" + k)(u + v) + {v') 2 + 2kv. 
Replacing u and v one obtains 

(v') 2 = S(x)(u + v) 2 - 2(u" + k)(u + v) + {u') 2 + 2ku. 



Subtracting one equation from the other one and using (|12|) one obtains 5 = 8 = 
const . Summing and dividing by ii + u give u" + v" = S(u + v) — k. Separation of 
variables yields (u 1 ) 2 — Su 2 + 2au + e, (v 1 ) 2 = 5v 2 + 2(3v + e, where a + (3 = —k, 
and substitution into ( |l2| ) proves e = e. The last of the equations ( |ll| ) is obtained 
in virtue of symmetry of the coordinate axes. I 

It is clear that solutions of the equations ( |ll| ) must satisfy also some addi- 
tional relations. However their analysis is not in principle difficult, and direct 
examination of all solutions brings to the following list. 

Theorem 2 The equations (j^j integrable in the sense of Definition 1, are ex- 
hausted, up to changes q rn . n = c(q m ^ n — mxo — m/o) o-nd permutations of coor- 
dinate axes, by the following sets of the functions f,g,h. In formulae (A), (B), 
( C) parameters are constrained by relation A + /i + v = 0, and in (I) by relation 
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Legendre transformations ftnfc together the equations corresponding to 

solutions (B) and (C), (D) and (E), (F) and (G), while equations corresponding 
to solutions (A), (H) and (I) are self-dual. I 

Notice, that cases (A) and (B) are connected by point transformation q = 
exp(2g). It is explained by fact that the Lagrangian ^(/ilogx + ^logy + Alogz) 
of the equation (||) , (A) is invariant under the dilations q — > Cq as well as under 
the shifts q — > q + C. From the other hand, inversions q — > q/(l — Cq) preserves 
Lagrangian as well but the change q = l/q which maps this group into the shift 
group does not bring to new equation. 

Apparently the cases (E), (F), (G), (H) and (A), (B) for A = appeared in 
the papers pi |L m, 0] by Suris for the first time. The rest cases are new, up to 
the author knowledge. 



3 Higher symmetries 

Study of the properties of the presented equations goes out of the scope of this 
paper. Of course, one has to prove that the Definition 1 brings to equations 
which are integrable in more habitual sense, that is satisfying the Painleve test or 
possessing the higher symmetries and conservation laws. Now we restrict ourself 
by presenting of the zero curvature representation and some higher symmetries 
for the most simple equation (||), (A), that is 

u a v v A A , 

= 0, (13) 



9i,o- 9 9-1,0 - 9 9o,i-9 90,-1 - 9 9i,i - 9 9-1.-1 

where A + /i + v = 0. Notice, that all formulae in this section remain valid also 
for vanishing A when the equation (^) degenerates into equation of the form 
(§■ 

(Continuous) symmetry of the discrete equation E = is a vector field 
qt = $ on the lattice which preserves this equation: dt(E) = 0\e=o- At first let 
us consider, for completeness, the classic symmetries q t = 1, qt = 9, qt = Q 2 
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corresponding to the group of linear-fractional transformations. Since these 
symmetries are variational, hence in virtue of the Noether theorem multiplying 
of ( p"3| ) by their right hand sides yields some conservation laws. Indeed, the 
shift symmetry generates the momentum conservation law (|^), and dilations 
and inversions give respectively the conservation laws 

(T m -1) m +(T W -1) Vq +(T m T n -l)- M -'■ 



9-9-1,0 q-qo,-i 9-9-1,-1 

(T m - D^tlSL + {Tn _ x) m^ + {TmTn _ 1} jfew = o. 
9-9-i,o 9-90,-1 9-9-1,-1 

The structure of the higher symmetries of the equation is rather orig- 
inal, as one can infer out of the simplest representatives, and it would be very 
interesting to obtain the complete description of the symmetry algebra. 

Theorem 3 The following formulae define the symmetries of the equation ([J 
1 fi v A 1/i v A 



9? 9-i,o-9 9o,i- 9 9-1,-1 - 9 9r, 9i,o- 9 90,-1 - 9 9-1,-1 -9 

1 11 v A 
— = - + + . 

9C 9i,o - 9 9o,i - 9 9i,i - 9 

Vector fields d^,d v ,d^ commute in virtue of the equation fil^j). I 

The proof can be obtained by straightforward, although tedious computa- 
tion. It is sufficient to check the determining equation dt(E) = 0|_e = o only for 
one symmetry, because two others are obtained by cyclic permutation of the 
axes x, y, z. 

The zero curvature representation for the equation ( [l3|) is defined as compati- 
bility condition for the auxiliary linear problems $1,0 = A^> and $x,i = ■B^'i.o- 
In order to retain the equal rights of the directions on the lattice, let us consider 
also the equation W = C^>\ t \ and obtain 

BA = A , 1 B_ lt0 , AC = C li0 A 1>1 , CB = B^^C^. (14) 

Consistency with the problem of the form ^> t — corresponding to the con- 
tinuous symmetry, brings to the equations 

At = Ux,oA — AU, B t = Ui, 1 B-BU 1 ,o, C t = UC - CU Xt i. (15) 



Matrices A, B and C are of the form 

A = {k - n)I - vP{q, ©,-1) - XP(q, 9-i,-i), 

B = kl - iiP{q, 91,0) - AP(?, 91.1), 

C = (k + v)I - fiP(q, g_ 1>0 ) - vP(q, 90,1), 
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where k is spectral parameter, / is unit matrix and projector P is given by 
formula 

1 / -v -uv 

P{u,v) = . 

u—v\l u 

(the coefficients on / can be chosen more symmetric by shift k — > k + (fi — v)/3).) 
On the picture above each of the matrices A, B, C depends on the vertices 
of the triangle in which it is placed (the center of the hexagon corresponds 
to the variable q). Notice also that successive application of the operators 
A, B, C maps the wave function into itself up to the scalar factor, since CBA = 
k(k — fi)(k+ v)I. The check of the equivalence of fll3| ) and equations Jl^ ) is not 
difficult when using the properties 



P(u,v)P(p,q) = 



(u - q)(p 



-P{p,v), P(u, v) + P(v, u) = I. 



{u-v){p-q) 

The representations JT^) for the symmetry dt = Oq are given by the matrix 
U = fc _1 (C7 — |7), where U is the projector 

1 ( a 



U = 



+ 1-1,00 



1,0" 

1,0/3 



a = Xqq-i -l + W-i,-iQo,-i + v Qo,-iQ> P = Mo-i + M + "3-1,-1- 

It should be stressed that the symmetries presented in Theorem 3 are not 
integrable equations by themselves and it makes sense to consider them only 
together with the constraint (|l3|). In other words, they are 1+1 rather than 1+2- 
dimensional equations. Indeed, denoting u m — q m>n , v rn = q m . n _i and using 
equation (13) one can rewrite the symmetries d v , as the one-dimensional 
lattices 

I u v A 1 m v A . 

(16) 



Ul 



V-l 



v-i 



Ul 
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it,£ it_i — u v — u V—i — u ty V\ — v u — v ui — v 

(in order to write down in analogous form one should choose dynamical 
variables along the axis y or z.) 

Remarkable fact is that these lattices are equivalent to one of the relativistic 
lattices (0), namely to the lattice, corresponding to the isotropic Heisenberg 
model. 



Theorem 4 Elimination of the variables v m from the equations ( \lq ) brings to 
the lattice 



* (Ml — U)- [U — 1t_lJ z U\ — VL u — U_l 



Vector fields defined by equations and ^ly commute and variables u = u r , 
v = v. m satisfy the system (now we assume that /i ^ 0) 

u v< = — + — '- - -^Wl - AXu- — '— ^ 

u — v [i /i U (u — v) z 



v — u fi /u y (u — v) 



in virtue of these equations. I 

So, the relation between the continuous equations (|^) and their discrete 
analogues (|^) is more profound than it seemed at first sight. Remind, that 
the relation between the lattices (||) and equations of the form (f|) has been 
already established in paper Q . More precisely, it was demonstrated there that 
the composition of the Legendre transformations TZ 1 T + for relativistic Toda 
lattices is equivalent to the shift in Toda lattices (|l|) (however, the analogy with 
discrete equations is slightly failed at this point, since the composition TZ 1 T + of 
the transformations @ cannot be rewritten as equation of the form (||)) and the 
nonlinear superposition principle for the lattices of these two types brings to the 
equations (^). The example considered above demonstrates that the equations 
(^|) can be interpreted as nonlinear superposition principle as well, but for the 
pair of relativistic lattices. 

In conclusion of this section notice that, although the equation (|l^) was 
not subjected to the Painleve test, but the simplest examples of 2-periodic 
reductions (several variants are possible) demonstrate that it is Painleve integr- 
able even at A + /i + v ^ 0. Of course in this case the above zero curvature 
representations and symmetries are lacking. Possibly, this situation is analogous 
to the case of difference Hirota-Miwa and KdV equations which were verified 
by the Painleve test in the paper M. 
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4 Concluding remarks 



The author's main purpose was to demonstrate that condition of invariance 
with respect to the Legendre transformations allows effectively isolate the class 
of integrable equations. The obtained examples of difference equations at the 
triangular lattice provide the discrete analogues of relativistic Toda lattices and 
prompt expectation that the method of Legendre transformations suggested in 
is not just a trick and its scope is rather wide. The open problem is further 
generalizations on the new classes of equations, in particular multidimensional 
ones. 

As in continuous case the weakest point of the method is the conjecture 
about the shift invariance of the Lagrangian. The study of the relativistic Toda 
lattices [p| implies that the difference analog of the Landau-Lifshitz equation 
must exist which is characterized by absence of the classic symmetries. 

^From the other hand, it is not difficult to find some multifield generaliza- 
tions, although their classification is far from completeness, as well as in con- 
tinuous case. Now we present only the vector analog of the discrete Heisenberg 
equation (|l3|): 

/j,(T m - 1) 7-|2 + v(T n - 1)— = (/i + u)(T m T n - 1) r-i2 
Fr \V\ \ z \ 

where q 6 R d and \q\ denotes Euclidean norm. Legendre transformation T + is 
of the form 

v Vos i . ^ x + y 0>1 v x x + y 0A 
X = v- {fi + u)- — ■ Y-ifl = -M|-|2 + (M + ^Ji — : (I 

wo.ir \x + yo,i\ Fr \x + yo,i\ 

and, as in the scalar case, maps equation into itself. Apparently this exam- 
ple admits generalizations for arbitrary Jordan triple systems with invertible 
elements. 
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